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Abstract 

We prove existence of weak solutions for a diffuse interface model for 
the flow of two viscous incompressible Newtonian fluids in a bounded 
domain by allowing for a degenerate mobility. The model has been 
developed by Abels, Garcke and Griin for fluids with different den- 
sities and leads to a solenoidal velocity field. It is given by a non- 
homogeneous Navier-Stokes system with a modified convective term 
coupled to a Cahn-Hilliard system, such that an energy estimate is ful- 
filled which follows from the fact that the model is thermodynamically 
consistent. 
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1 Introduction 

Classically the interface between two immiscible, viscous fluids has been 
modelled in the context of sharp interface approaches, see e.g. |Mue85j . But 
in the context of sharp interface models it is difficult to describe topological 
changes, as e.g. pinch off and situations where different interfaces or different 
parts of an interface connect. In the last 20 years phase field approaches have 
been a promising new approach to model interfacial evolution in situations 
where interfacial energy effects are important, see e.g. [CheC)2]. In phase 
field approaches a phase field or order parameter is introduced which rapidly 



*Fakultat fiir Mathematik, Universitat Regensburg, 93040 Regensburg, Germany, e- 
mail: helmut.abels@mathematik.uni-regensburg.de 

^Fakultat fiir Mathematik, Universitat Regensburg, 93040 Regensburg, Germany, e- 
mail: daniel.depner@mathematik.uni-regensburg.de 

"^Fakultat fiir Mathematik, Universitat Regensburg, 93040 Regensburg, Germany, e- 
mail: harald.garcke@mathematik.uni-regensburg.de 



1 



changes its value in the interfacial region and attains two prescribed values 
away from the interface. 

For two-phase flow of immiscible, viscous fluids a phase-field approach 
first has been introduced by Hohenberg and Halperin [HH77j, the so-called 
"Model H". In their work the Cahn-Hilliard equation was coupled to the 
Navier-Stokes system in such a way that capillary forces on the interface are 
modelled with the help of the phase field. The approach of Hohenberg and 
Halperin [HH77] was restricted to the case where the densities of the two 
fluids are the same or at least are very close ("matched densities"). It has 
been later shown by Gurtin, Polignone, Vihals [GPV96j that the model can 
be derived in the context of rational thermodynamics. In particular global 
and local energy inequalities are true. These global energy estimates can be 
used to derive a priori estimates and this has been used by Boyer |Boy99| 
and by Abels |Abe09b| for proofs of existence results. 

Often the densities in two phase flow are quite different. Therefore, there 
have been several attempts to derive phase field models for two phase flow 
with non- matched densities. Lowengrub and Truskinovsky |LT98| derived 
a first thermodynamically consistent phase field model for the case of dif- 
ferent densities. The model of Lowengrub and Truskinovsky is based on 
a barycentric velocity and hence the overall velocity field turns out to be 
not divergence free in general. In addition, the pressure enters the Cahn- 
Hilliard equation and as a result the coupling between the Cahn-Hilliard 
equation and the Navier-Stokes equations is quite strong. This and the fact 
that the velocity field is not divergence free makes numerical and analytical 
approaches quite difficult. To the authors knowledge there have been so far 
no numerical simulations for the full Lowengrub- Truskinovsky model. With 
respect to analytical results we refer to the works of Abels |Abe09a[ IAbel2] 
for existence results. 

In a paper by Ding, Spelt and Shu |DSS07j a generalization of Model 
H for non-matched densities and a divergence free velocity field has been 
derived. However it is not known whether this model is thermodynami- 
cally consistent. A first phase field model for non-matched densities and 
a divergence free velocity field which in addition fulfills local and hence 
global free energy inequalities has been derived by Abels, Garcke and Griin 
|AGG12j . The model in |AGG12) is given by the following system of Navier- 
Stokes/Cahn-Hilliard equations: 

dt{p{ip)v) + div(v (g) {p{ip)v + J)) - div(2T/((/?)Dv) + Vp 

= — div{a{(p)Vip (8) Vtp) in Qt, 



div V 







in Qt, 
in Qt, 



dt(p + V • Vip 



div {m{ip)ViJ,) 




div {a{ip)Vip) 



in Qt, 
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where J = -^^m((/?)V/i, Qt = ^ x (0,T) for < T < oo, and ft C M'', 
d = 2, 3, is a sufficiently smooth bounded domain. We close the system with 
the boundary and initial conditions 

v|9n = ondQxiO,T), 
dnV>\dn = dnfJ-lan = on dil x (0, T), 

(v,v9) |t=o = (vo,<^o) in J^, 

where = n ■ and n denotes the exterior normal at dQ. Here v is the 
volume averaged velocity, p = p{ip) is the density of the mixture of the two 
fluids, (/? is the difference of the volume fractions of the two fluids and we 
assume a constitutive relation between p and the order parameter ip given 
by p{(p) = ^{pi + P2) + ^{P2 — Pi)^, see [ADG12j for details. In addition, 
p is the pressure, p is the chemical potential associated to ip and pi, p2 
are the specific constant mass densities of the unmixed fluids. Moreover, 
Dv = i(Vv + Vv^), r]{ip) > is a viscosity coefficient, and m{ip) > is a 
degenerate mobility coefficient. Furthermore, ^{p) is the homogeneous free 
energy density for the mixture and the (total) free energy of the system is 
given by 

Eireei^) = ^ (^(</.) + a((/.)^) dx 
for some positive coefficient a{ip). The kinetic energy is given by E]^^^{ip, v) = 

I 1 2 

In p{y^)'2 ^^^'^ total energy as the sum of the kinetic and free energy 
is 

= ^ dx + J^ + aM^) dx. ^^'^^ 

In addition there have been further modelling attempts for two phase flow 
with different densities. We refer to Boyer |Boy02| and the recent work of 
Aki et al. |ADGK12] . We remark that for the model of Boyer no energy 
inequalities are known and the model of Aki et al. does not lead to velocity 
fields which are divergence free. 

In |ADG12j an existence result for the above Navier-Stokes/Cahn-Hil- 
liard model has been shown in the case of a non-degenerate mobility m{ip). 
As is discussed in [AGG12j the case with non-degenerate mobility can lead 
to Ostwald ripening effects, i.e., in particular larger drops can grow to the 
expense of smaller ones. In many applications this is not reasonable and as 
pointed out in |AGG12j degenerate mobilities avoid Ostwald ripening and 
hence the case of degenerate mobilities is very important in applications. In 
what follows we assume that m{Lp) = 1 — 99^ for < 1 and extend this by 
zero to all of M. In this way we do not allow for diffusion through the bulk. 
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i.e., the region where ip = 1 resp. if = —1, but only in the interfacial region, 
where \ip\ < 1. The degenerate mobihty leads to the physically reasonable 
bound Iv?! < 1 for the order parameter ip, which is the difference of volume 
fractions and therefore we can consider in this work a smooth homogeneous 
free energy density ^ in contrast to the previous work |ADG12) . 

For the Cahn-Hilliard equations without the coupling to the Navier- 
Stokes equations Elliott and Garcke [EG96] considered the case of a de- 
generate mobility, see also Griin |Gru95| . We will use a suitable testing 
procedure from the work |EG96j to get a bound for the second derivatives 
of a function of (p in the energy estimates of Lemma 13.71 We point out that 
our result is also new for the case of model H with degenerate mobility, 
i.e., pi = p2, which implies J = in the above Navier-Stokes/Cahn-Hilliard 
system. 

The structure of the article is as follows: In Section 2 we summarize 
some notation and preliminary results. Then, in Section 3, we reformulate 
the Navier-Stokes / Cahn-Hilliard system suitably, define weak solutions and 
state our main result on existence of weak solutions. For the proof of the 
existence theorem in Subsections 3.2 and 3.3 we approximate the equations 
by a problem with positive mobility and singular homogeneous free en- 
ergy density For the solution (ve,(/?£, J^) of the approximation (with 
Je = —mi;{{pi;)'V^s) we derive suitable energy estimates to get weak limits. 
Then we extend the weak convergences to strong ones by using methods 
similar to the previous work of the authors |ADG12j . careful estimates of 
the additional singular free energy density and by an additional subtle ar- 
gument with the help of time differences and a theorem of Simon jSim87j . 
We remark that this last point would be easier in the case of a constant 
coefficient a{(p) in the free energy. Finally we can pass to the limit e — t- in 
the equations for the weak solutions (ve,(/9e, J^) and recover the identities 
for the weak solution of the main problem. 

2 Preliminaries and Notation 

We denote a (g) 6 = {aibj)fj^^ for a,6 e M'^ and ylsym = ^(^ + -4"^) for a 
matrix yl G M"'^^. If X is a Banach space and X' is its dual, then 

{f,9)^{f,9)x',x = f{g), fex',gex, 

denotes the duality product. We write X y if X is compactly em- 

bedded into Y. Moreover, if is a Hilbert space, (• , ■)h denotes its inner 
product. Moreover, we use the abbreviation (. , ■)m = {■ , ■)l^{m)- 

Function spaces: If M C M'^ is measurable, L'^{M), 1 < q < oo, denotes 
the usual Lebesgue-space and ||.||g its norm. Moreover, L^{M;X) denotes 
the set of all strongly measurable g-integrable functions if g G [l,oo) and 
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essentially bounded strongly measurable functions, if q = oo, where X is a 
Banach space. 

Recall that, if X is a Banach space with the Radon-Nikodym property, 
then 

L'?(M; Xy = L'^'{M; X') for every 1 < g < oo 

by means of the duality product {f,g) = Jj^^{f{x), g{x))x' ,xdx for / € 
L9'(M; X'), g G L'?(M; X). If X is reflexive or X' is separable, then X has 
the Radon-Nikodym property, cf. Diestel and Uhl [DU77] . 

Moreover, we recall the Lemma of Aubin-Lions: If Xq Xi ^ X2 

are Banach spaces, 1 < p < 00, I < q < 00, and / C M is a bounded interval, 
then 

E LP{I;Xo) : ^ G L''(/;X2)} LP{I;X^). (2.1) 

See J.-L. Lions |Lio69j for the case q > I and Simon |Sim87j or Roubi- 
cek |R,ou9nj for q = 1. 

Let 17 C M'^ be a domain. Then W^{n), k £ Nq, I < q < 00, de- 
notes the usual L'^-Sobolev space, WI^q{Q) the closure of C^(il) in Wg{n), 
W-''{n) = {W^, Q{n)y, and W-^{n) = {W^,{n)y. We also use the abbrevi- 
ation H''{n) = w^{n). 

Given / S L^{Q), we denote by /n = py f^f{x)dx its mean value. 
Moreover, for m € M we set 

q^^in) := {/ G L'^in) ■.fn = m}, 1 < < cx). 
Then for / G L^(il) we observe that 

Pof :=f-h = f-^J'^f{x)dx 

is the orthogonal projection onto L^Q^(il). Furthermore, we define 

Then H^q^{^) is a Hilbert space due to Poincare's inequality. 

Spaces of solenoidal vector-fields: For a bounded domain C M"^ we 
denote by C^(r2) in the following the space of all divergence free vector 
fields in Cq°{Q)'^ and L'^{^) is its closure in the L'^-norm. The corresponding 
Helmholtz projection is denoted by P^, cf. e.g. Sohr [SohOlj . We note that 
Paf = f — Vp, where p S W2{^) H L^Q^(il) is the solution of the weak 
Neumann problem 

{Vp,Vip)n = {f,Vip) for all 99 G (2.2) 
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Spaces of continuous vector-fields: In the following let / = [0, T] 
with < T < oo or let / = [0, oo) if T = oo and let X be a Banach 
space. Then BC{I;X) is the Banach space of all bounded and continu- 
ous f : I ^ X equipped with the supremum norm and BUC{I;X) is the 
subspace of all bounded and uniformly continuous functions. Moreover, we 
define BCw{I] X) as the topological vector space of all bounded and weakly 
continuous functions / : I ^ X. By C^(0, T; X) we denote the vector space 
of all smooth functions /: (0, T) — )• X with supp/ CC (0, T). We say that 
/ G H^pHO, T- X) for 1 < p < oo, if and only if /, f G ^^(0, T; X), where f 
denotes the vector- valued distributional derivative of /. Finally, we note: 

Lemma 2.1. Let X, Y be two Banach spaces such that Y X and X' ^ Y' 
densely. Then L°°{I;Y) r\ BUC{r,X) ^ BCy,{I;Y). 

For a proof, see e.g. Abels |Abe09a| . 



3 Existence of Weak Solutions 

In this section we prove an existence result for the Navier-Stokes/Cahn- 
Hilliard system from the introduction for a situation with degenerate mo- 
bility. Since in this case we will not have a control of the gradient of 
the chemical potential, we reformulate the equations by introducing a flux 
J = —m{ip)V consisting of the product of the mobility and the gradient of 
the chemical potential. In this way, the complete system is given by: 

5t(pv) -I- div(pv (g) v) — div(2?7(99)L'v) + Vp 

+ div(v o /3J) = - div(a((^)V(/9 in Qt, (3.1a) 

divv = in Qt, (3.1b) 

SfV? -|- v • = — div J in Qt, (3.1c) 

J = -m(^)V (^vI/V)+aV)^ 

-div{a{(f)V(f)^ in Qt, (3.1d) 

v|an = on St, (3.1e) 

dny^ldn = (J • n)\dQ = on St, (3.1f) 

(v,(/?) \t=o = (vo,93o) in ^, (3.1g) 

where we set /? = and J = —m{ip)'V fi as indicated above. The 

constitutive relation between density and phase field is given by p{if) = 
^(Pi + P2) + ^{P2 — Pi)^ as derived in Abels, Garcke and Griin |AGG12| . 
where Pi > are the specific constant mass densities of the unmixed fluids 
and ip is the difference of the volume fractions of the fluids. By introducing 
J, we omitted the chemical potential p, in our equations and we search from 
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now on for unknowns (v,(/?, J). In the above formulation and in the fol- 
lowing, we use the abbreviations for space-time cylinders Q(s,t) = ^ x {s,t) 
and Qt = Q{o,t) and analogously for the boundary = x (s, t) and 
St = 'S'(o,t)- Equation (j3.1ep is the no-slip boundary condition for viscous 
fluids, (J • n)|an = resulting from dnfJ-\dn = means that there is no mass 
flux of the components through the boundary, and dn(p\dn = describes 
a contact angle of 7r/2 of the diffused interface and the boundary of the 
domain. 



3.1 Assumptions and Existence Theorem for Weak Solutions 

In the following we summarize the assumptions needed to formulate the 
notion of a weak solution of (j3.1ap -(3.1g) and an existence result. 

Assumption 3.1. We assume that Q, C M'^, d = 2,3, is a bounded domain 
with smooth boundary and additionally we impose the following conditions. 

(i) We assume a, ^' G C^M), r/ G C°(IR) and < cq < a{s),ri{s) < K for 
given constants cq,K > 0. 

(ii) For the mobility m we assume that 

We remark that other mobilities which degenerate linearly at s = ±1 
are possible. The choice (j3.2|) typically appears in applications, see Cahn 
and Taylor |CT94] and Hilliard [Hil70) . Other degeneracies can be handled 
as well but some would need additional assumptions, see Elliott and Garcke 
[EG96]. 

We reformulate the model suitably due to the positive coefficient a((p) in 
the free energy, so that we can replace the two terms with a{ip) in equation 
()3.1dp by a single one. To this end, we introduce the function A(s) := 
Jo y(i{j)dT. Then A'{s) = ya{s) and 



-M^)AA{^) = a'iip) ^ - div (a((^) Vip) 

resulting from a straightforward calculation. By reparametrizing the poten- 
tial through $ : M M, $(r) := ^'(yl"^(r)) we see ^''(s) = a{s)^' {A{s)) 
and therefore we can replace line ()3.1dp with the following one: 

J = -m^V (v/^ (^'(^M) - A^M)) . (3.3) 
We also rewrite the free energy with the help of A to 

^free(v^) = (^(^(V^)) + dx . 
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Remark 3.2. With the above notation and with the calculation 
— div(a((/9)V(/7 Vip) 

= - dw{a{ip)Vip)Vip - a{ip)V 
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we rewrite line (|3.1ap wzi/i a new pressure g = p + a{Lp)^-^Y~ into: 

dt{p\r) + div(pv (g) v) - div(2?7((^)i:>v) + + div(v (g) /3J) 

I (3-4) 

We remark that in contrast to the formulation in \ADG1^ we do not use 
the equation for the chemical potential here. 

Now we can define a weak solution of problem (j3.1ap - p.lgp . 

Definition 3.3. Let T £ (0,oo), vq G LI{U) and (po G H^{Q) with |(/7o| < 1 
almost everywhere in Q. If in addition Assumption \3.1\ holds, we call the 
triple (v, (f, J) with the properties 

V G BC^iiO, T];Ll{n)) D L\0, T; H^{nf) , 

(f G BC^{[0,T];H^{n))nL'^{0,T;H^{n)) with \ip\ < 1 a.e. in Qt , 
J G L^{ld,T;L^(Slf ) and 

(v,V3) \t=Q = (vo,V3o) 

a weak solution of ()3.1aP - p.lgP if the following conditions are satisfied: 
- (P^, dtip)Q^ + (div(pv v), + i2ij{^)Dv, DxP)^ 

- ((v /3J), Vt/.)Q^ = - (y^Ayl((^) V(/p, 

for all -0 G [C^(n X (0,T))]'^ with div^/^ = 0, 

-/ ipdtC,dxdt+ i ■Vip)Cdxdt = / J-VCdxdt (3.6) 

Jqt Jqt Jqt 

for all C G C^((0,r;Ci(lI)) and 



(3.5) 



J • ?7 (it 

(^ya(^ (^$'(^((/p)) - A^(99)) ) div(m(v9)r7) 



(3.7) 



for all rj G L'^{0, T; H^{n)'^) n L°°{QtY which fulfill r]-n = on St- 
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Remark 3.4. The identity ()3.7p is a weak version of 



Our main result of this work is the fohowing existence theorem for weak 
solutions on an arbitrary time interval [0, T], where T > 0. 

Theorem 3.5. Let Assumption \3.1\ hold, vq E L'^{^) o^f^d £ H^{Q) 
with \{Pq\ < 1 almost everywhere in fi. Then there exists a weak solution 
(v,ip,J) of ()3.1ap - p.lgP in the sense of Definition \3.3[ Moreover for some 



J G L'^{Qt) it holds that J = y^m{if)J and 

Etot{'p{t),v{t)) + 2^{ip)\Dw\^ dxdr + I \3\'^dxdT 

•^Q{3,t) •^Q{s,t) (3-8) 

<Etot{^{s),v{s)) 

for all t G [s,T) and almost all s £ [0,T) including s = 0. The total energy 
Eiot is the sum of the kinetic and the free energy, cf. (jl.ip . In particular, 
J = a.e. on the set {\ip\ = 1}. 

The proof of the theorem will be done in the next two subsections. But 
first of all we consider a special case which can then be excluded in the 
following proof. Due to \(Pq\ < 1 a.e. in Q we note that -f^ipodx S [—1, 1]. 
In the situation where -j-^ipodx = 1 we can then conclude that tpQ = \ 
a.e. in Q. and can give the solution at once. In fact, here we set = 1, 
J = and let v be the weak solution of the incompressible Navier-Stokes 
equations without coupling to the Cahn-Hilliard equation, where p and rj are 
constants. The situation where -j-^^p^dx = —1 can be handled analogously. 

With this observation we can assume in the following that 

f (fodx e (-1,1) , 
Jn 

which will be needed for the reference to the previous existence result of the 
authors |ADG12| and for the proof of Lemma \37l\ (Hi). 

3.2 Approximation and Energy Estimates 



In the following we substitute problem (|3.1ap -( 3.1gP by an approximation 



with positive mobility and a singular homogeneous free energy density, which 
can be solved with the result from the authors in jADG12j . For the weak 
solutions of the approximation we then derive energy estimates. 

First we approximate the degenerate mobility m by a strictly positive 

me as 

m(— 1 + e) for s<— 1 + e, 
m^(s) := ( m{s) for |s| < 1 — e , 
m(l — e) for s > 1 — e . 
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In addition we use a singular homogeneous free energy density given by 

^e{s) ■■= ^(s) + e^inls) 1 where 
^ln(s) := (1 + s) ln(l + s) + (1 - s) ln(l - s) . 

Then G C([-l, 1]) Pi C^((-l, 1)) fulfiUs the assumptions on the homo- 
geneous free energy as in Abels, Depner and Garcke |ADGl2j . which were 
given by 

lim "^'(s) = ±oo , ^'l'(s) > K for some k € M and lim ! = +oo. 

s^±l '^^ ' e\ ^ - s^±l ^^s) 

To deal with the positive coefficient a{(p), we set similarly as above ^in('') '■= 
and ^e{r) := ^'^(^-^(r)) for r G [a,b] := 
Now we replace m by and ^' by and consider the following ap- 
proximate problem, this time for unknowns {v,ip,fi): 

dt{pv) + div (pv (g) v) - div {2ri{ip)Dw) + V5 

-hdiv(v(g/3me((/?)V/i) = -ya(^AA(93)Vv7 in Qt, (3.9a) 
divv = in Qt, (3.9b) 

(?f(/7 -|- V • V(/? = div(?Tie((/9)V;u) in Qt, (3.9c) 



/i = (^^(^(9^)) - AA((^)j in Qt, (3.9d) 

v|an = on St, (3.9e) 

dnf\dn = dnnlan = on St, (3.9f) 

{\r,ip) \t=o = (vo,(/?o) in ^- (3.9g) 

Prom [ADGI2) we get the existence of a weak solution (v^, ^u^) with the 
properties 

^r,eBCU[0,T];Ll{n))nL\0,T■,H^{ny), 

iPseBCU[0,nH\n))nL\0,T;H\Q)), ^',{ip,) e L\0,T; L\n)) , 
lie G L'^{0,T;H^{n)) and 
(v£,V'£)lt=o = (vo,y'o) 
in the following sense: 



(pgVg, 9ti/')Qj, + (div(p£V£ Ve), -0)^^ {2r]{ipe)D^re, Dtl)) 



T 

T 

d 



(3.10) 



for ah -0 G [C^(f^ X (0,T))]'' with div-0 = 0, 

- iv^e, dtOg^ + (v, • V^s, Oq^ = - (m,((/p,)V/x„ VOq^ (3.11) 
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for all C G C^{{0,T);C\n)) and 

fis = V o.{(pe) (^^'^{A{ipi;)) — AA{ipi;)^ almost everywhere in Qt- (3.12) 
Moreover, 

EtotiMt),Mt))+ [ 2r]{ipe)\DVe\'^dxdT 

y^'" (3.13) 



+ 



for all t £ [s,r) and almost all s S [0, T) has to hold (including s = 0). 

Herein is given as = \{pi + P2) + \{P2 — Pij^e- Note that due to 
the singular homogeneous potential ^'^ we have \tpe\ < 1 almost everywhere. 

Remark 3.6. Note that equation ()3.10p can be rewritten with the help of 
the identity 



This can he seen by testing p.l2p with Vfe ■ V' and noting that is diver- 
gence free. 

For the weak solution (v^, //g) we get the following energy estimates: 

Lemma 3.7. For a weak solution (v^, (/j^, /Xg) of problem ()3.9ap -(3.9g) we 

have the following energy estimates: 

(i) sup / (p^(t)^^^ + hVipe{t)\^ + ^eiMt)))dx 

o<t<TJn \ I ^ J 

+ / 2ri{ipe)\Dwe\'^ dxdt+ / me{^e)\^ P-e? dx dt < C, 
Jqt JQt 

(a) sup / Gs{fs{t))dx+ / \AA{ips)\^dxdt < C, 

0<t<TJn JQt 
{Hi) £^ I \^'i^{ip,)\^dxdt < C, 

Jqt 

Here Gs is a non-negative function defined by Ge(0) = 0^(0) = and 
G':{s) = ^^y^) for se [-1,1]. 

Proof, ad (i): This follows directly from the estimate ()3.13p derived in the 
work of Abels, Depner and Garcke |ADG12j . We just note that for the esti- 
mate of V(/9e we use VA{ifir) = y^a(Jp^Vipe and the fact that a is bounded 
from below by a positive constant due to Assumption 13. 1[ 
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ad (n): From line ([XTT]) we get that dtife G L^{0,T; (H^in))'), since 
G LP'{Qt) and v • Vip = div(v(/?) with v(/? G L'^{Qt)- Then we derive 
for a function ^ G L'^{0,T; H'^{Q,)) the weak formulation 

/ C)dT + • V(/3£:C fi^a; dr 

JO JQt 

= -/ me{(fs)^fJ'e -^CdxdT (3.14) 
= ^ - AAi^e)) div(m,(99e)VC) dr, 

where we additionally used ()3.12p to express fig. Now we set as test function 
( = G',{ipe), where is defined by Ge(0) = G',{0) = and G'^{s) = 
^^^^^ A'{s) for s G [—1,1]. Note that Gg is a non-negative function, which 

can be seen from the representation Gs{s) = ^ j^^^j^A'{t) dr^ dr. With 
( = G'g{(pe) it holds that 



VC = G''{ipe)V(ps = {A{ipe)) and therefore 



div(m,(v9,)VC) = A(A(v9,)). 



Hence we derive 

ft 



[ {dtipe,G',{^e))dT + [ Ye-y^eG',{^e)dxdT 
JO JQt 

= (^'MiVe)) - AA(^,)) AA{^,) dx dr (3.15) 

■^'g{^e)AA{^e)dxdT - / ^Ja{<pe) \AA{^e)\^ dx dr. 
t JQt 



With this notation we deduce 
rt 



{dt(p£,G'g{ipe))dt = / Ge{<p{t))dx - / Ge{<PQ)dx and 
Jn Jn 

/ ^£ ■VipeG'g{'fe)dxdt = / We -V {Ge{<P£))dxdt 
JQt JQt 

= — I d\Y^^£Ge{(pe) dx dt = Q . 
JQt 

For the first term on the right side of (j3.15p we observe 

%{<Pe)AA{<fe)dx dr 

= / ^'{ipe)AA{ip,)dxdT + e / ^[^{ipe)AA{ip,)dxdT 
JQt JQt 
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< - / ^"i^e)^Ve ■ VA((/J^) dx dt 
jQt 

= - I ^"{^e)Va{Ve)\^Ve\'^dxdt. 
JQt 

Herein the estimate 

/ ^[^{^,)^A{^,)dxdT <{) 

JQt 

for the logarithmic part of the homogeneous free energy density is derived 
as follows. With an approximation of ipe by (/j" = for < a < 1 we 
have that < a < 1 and therefore 

/ ^;j(/p-)A^((^°)dxdr = - / vI/^;((^^)V(^^Vyl((^^)dx(^r <0, 
JQt JQt 

where we used integration by parts. To pass to the limit for a 1 in the left 
side we observe that — t- in L^(0, T; //^(fi)). Hence together with the 
bound < l^}n(¥'e)l can use Lebesgue's dominated convergence 

theorem to conclude 

/ ^;j(^°)AA(^-)dxdr ^ / ^[^[pe)^A{p,)dxdT for a I. 
JQt JQt 

With the bound from below a(s) > cq > from Assumption 13 . 1 1 we derived 
therefore 

f Ge{p{t))dx+ f \^A{ipe)\^ dxdT 

<C( [ Ge{po)dx+ [ ^"{pe)Va{^e)\^Vs\^dxdT) . 

\Jn JQt J 

Now we use m^{T) > m{T) to observe the inequality 




Due to the special choice of the degenerate mobility m in (|3.2p we conclude 
that G can be extended continuously to the closed interval [—1,1] and that 
therefore the integral G{ipo) dx and in particular the integral Ge(v9o) dx 
is bounded. 

Moreover, since ^"{s) is bounded in |s| < 1 and since we estimated 
Jq |V<y9e(t)p dx in (i), we proved (ii). 
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ad (iii): To show this estimate we will argue similarly as in the time- 
discrete situation of Lemma 4.2 in Abels, Depner and Garcke |ADG12) . 
We multiply equation p.l2p with PoV's) integrate over Q and get almost 
everywhere in t the identity 

/ HePo^£dx= / '^'{ipe)Po^edx + £ / '$[^{ipe)Po(Pe dx 

Jfi Jo. Jq ^2 -j^g^ 

y/a{(pe)AA{ipe)Po(Ps dx. 



n 



By using in identity (j3.11|) a test function which depends only on time t and 
not on X £ we derive the fact that {(Pe)n = {^o)n and by assumption this 
number lies in (— l+a, 1—a) for a small a > 0. In addition with the property 
lims_j,-|-i ^'j^(s) = ±oo we can show the inequality ^'j^(s)(s — {(po)n) > 
Cq,|^'Jj^(s)| — Cq, in three steps in the intervals [— 1, — 1 + [—1 + ^,1 — ^] 
and [1 — ^,1] successively. Altogether this leads to the following estimate: 

e [ \^[Jve)\ dx<c(e [ ■i'[Jipe)Po^e dx + l) . (3.17) 
Jn \ Jn J 

We observe the fact that /igPoV'e dx = j^{PQ^e)'^£ dx and due to integra- 
tion by parts 

- / a/ a{Lpe)AA{Lpe)PQiPe dx 

Jn 

= / V a{ip^)V A{Lpe) ■ Vife dx + / -a{Lpey'^\/ipe-VA{ipe)Poipedx 

Jn Jn ^ 

= / a{(ps)\VLpe\'^ dx + / -Po(pe\VLpe\'^ dx . 

Jn Jn ^ 

Combining estimate (j3.17p with identity (|3.16p we are led to 

el mj^e)\dx<c( [ \{Pofle)^e\dx+ [ i^,)Po^e\ dx 

Jn \Jn Jn 

+ Jjy/ai^e)AA{^e)Po^e\ dx + 1^ 

<C{\\Pofie\\LHn) + \N^e\\mn) + '^) 
<C{\\Vfie\\mn) + l)- 

In the last two lines we have used in particular the facts that ipe is bounded 
between —1 and 1, that ^' is continuous, the energy estimate from (n) for 
supo<t<r ||V99e||/^2(Q) and the Poincare inequality for functions with mean 
value zero. 
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With the last inequahty we can estimate the integral of fis by simply 
integrating identity ()3.12p over Q: 





<-l 


Jn 


Jn 



\^'{ipe)\dx + e / mjyDe)\dx + 
Jn 



'\/a{ipe)AA{ips) dx 



n 



<C(||V^e||L2{Q) + l), 

where we used similarly as above integration by parts for the integral over 
a{^pe)A.A{ipi.). By the splitting of fis into /i^ = Po/"e + {fJ's)n we arrive at 

Mh^n) < C (\\Vfie\\l2^n) + l) • 
Then, again from identity ()3.12p . we derive 

and together with the last estimates and an additional integration over time 
t this leads to 

£'ll*'m(^.)lli^(Q,)<^^(l|V^.lli^(Q,) + l)- 

Note that we used the bound ||A^((^£)||2,2(-Qy) < C from (ii). Furthermore, 
due to the bounds in (i), we see £||V//£||^2(Qj,-) < C since m^^s) > e for 
|s| < 1 and therefore we arrive at 

ad (iv): This follows directly from (i). □ 

3.3 Passing to the limit in the Approximation 

In this subsection we use the energy estimates to get weak limits for the se- 
quences {vs,iPe,Je), where = ^/mJJpsjJe (= -ms{ips)y f^e)- With some 
subtle arguments we extend the weak convergences to strong ones, so that 
we are able to pass to the limit for e — )• in the equations (j3.10p - (j3.12p to 
recover the identities (|3.5p - (|3.7p in the definition of the weak solution for 
the main problem (j3.1ap-([3.1gp. 



Using the energy estimates in Lemma 13.71 we can pass to a subsequence 
to get 

^ V in L'^{0,T;H'^{ny), 
^e^^ in L\0,T;H\n)), 
Je^J in L'^{0,T-L^{n)'^) and 
Je^J in L'^{0,T;L^{n f) 
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for V G L^{0,T; H\ny)nL°°{0,T; Ll{n)), ip G L°°(0,T; /^^(ri)) and J, J G 
L^(0, T; L-^(r2)'^). Here and in the following all limits are meant to be for 
suitable subsequences — )■ for A; — )• oo. 

With the notation = —mir{(pi;)V fie the weak solution of problem 
()3.9ap - p.9gP fulfills the following equations: 

(3.18) 



- ((v, ®/3J,),V'0)q^ = - (^^/a{ips)AAi^s)'^Vs,'4^ 
for all G [C^{^ X (0,r))]'^ with div-0 = 0, 

-/ (pedtCdxdt+ / ■V(pe)Cdxdt = / Js-VCdxdt (3.19) 

Jqt Jqt Jqt 

for ah C e C^((0,r;Ci(lI)) and 
Je • ridx dt 



Q 



T 



T 



^'^(99e) - '\/a{(ps)AA{ips)) div{me{(pe)ri) dx dt 



(3.20) 



for all T] G L2(0, T; H^{ny) n L°°{QtY with 77 • n = on 5t. For the last 
line we used that for functions with • n = on St it holds 

/ ie''ndxdt= I V p^ ■ m^{ip^)'q dx dt = — I jj^ dw{mi;{ipi;)r]) dx dt 
Jot Jot Jqt 

%{'Pe) - \/a{^e)^M'^e)] d\v{me{ipe)'n) dx dt . 



'T 



Now we want to pass to the limit e — )■ in the above equations to achieve 
finally the weak formulation (j3.5p - (j3.7p . 

For the convergence in identity ()3.18p we first note that 

dtipe is bounded in L^{Q,T]{H^{n))') and 
ipe is bounded in L°^(0, T; iJ^(17)) . 

Therefore we can deduce from the Lemma of Aubins-Lions ()2.ip the strong 
convergence 

tpe^ip in L^{id,T;L^{Q)) 

and — )• pointwise almost everywhere in Qt- 
From the bound of Aj4((/9e) in LP'{Qt) and from 

V A{lp^) ■ n = \J a{ipe)V^e ■ n = Q on St, 
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we get from elliptic regularity theory the bound 



ll^(¥'e)llL2(0,T;H2(n)) < C . 

This yields 

A{^,)^g in L\Q,T;H^{n)) 

at first for some g G i^(0, T; H'^{Q)), but then, due to the weak convergence 
S/ife Vip in L^(0, T; L^(r2)) and due to the pointwise almost everywhere 
convergence a{ipe) a{if) in Qt we can identify g with A{ip) to get 

^((/j,) ^ in L\0,T-H\n)). 

The next step is to strengthen the convergence of Vc/J^ in LP'{Qt)- To this 
end, we remark that by definition A is Lipschitz-continuous with 



\A{r)-A{s)\< 



\l a(r) dr 



< C\r- 



Furthermore from the bound of dt^^ in -^^^(0,T; (//^(ri))') we get with the 
notation + /i) for a shift in time 

+ /i) - ¥'£|lL2(o,'r_;j.(j^i(n))') < C/i , 

which leads to the estimate 

L2(0,T~/i;{_f/i{f^))') 

< C||99e(. + /l) - 92e|lL2(0,T-h;{Hi(!^))') 

< C/i — > as /i . 

Together with the bound of A[^^) in L2(0,T; we can use a theorem 

of Simon [Sim871 Th. 5] to conclude the strong convergence 

A{>p,)^A{^) in L\Q,T-H\^)). 



From V^(v9e) = ^Ja{}p£)Vipe we get then in particular the strong conver- 
gence 

V(^e^V99 in L2(0,r;L2(Jl)). 

In addition we want to use an argument of Abels, Depner and Garcke from 
|ADG12[ Sec. 5.1] which shows that due to the a priori estimate in Lemma 
13.71 and the structure of equation (|3.18|) we can deduce the strong conver- 
gence Vg — )• V in L'^(0, T; L^(O)'^). In few words we show with the help 
of some interpolation inequalities the bound of dtiPaiPe^e)) in the space 
L^W^iQ)') and together with the bound of Paipe'^e) in L'^{0,T; H^{nY) 
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this is enough to conclude with the Lemma of Aubin-Lions the strong con- 
vergence 

P,(p,v,)^P,(pv) in L^{Q,T-L^{^f). 

From this we can derive — )• v in L^(0, T; L^(i7)"'). For the details we refer 
to |ADG12l Sec. 5.1 and Appendix]. 

With the last convergences and the weak convergence Jg ^ J in LP'{Qt) 
we can pass to the limit e — t- in line ()3.18p to achieve (jS.Sp . 

The convergence in line (|3.19|) follows from the above weak limits of 
and Je in LP'{Qt) and the strong ones of and V(^£ in LP'{Qt)- 

Finally, the convergence in line (j3.2Up can be seen as follows: The left 
side converges due to the weak convergence of and for the right side we 
calculate 



^'^((/Je) - ^/ a{ipe)l^A{ipe)j d\Y{me{ipe)'n) dx dt 

'I^'{iPe)div{me{(pe)r])dxdt + £ / ^[^{ips) div{ms{(pe)'n) dx dt 

JQt 

- I ^/a(^AA{^ps)div{■me{(pe)r|)dxdt. (3.21) 
JQt 

The first and the third term can be treated similarly as in Elliott and Garcke 
|EG96] . For the convenience of the reader we give the details. 

First we observe the fact that m-g — t- m uniformly since for all s G R it 
holds: 

|me(s) — m{s)\ < m(l — e) — )• for e — )• . 

Hence we conclude with the pointwise convergence ip^ ^ <p a.e. in Qj- that 

ms{ipe) — > m{if) a.e. in Qt ■ 

In addition with the convergences "^'{ipe) — )• "if'i^p), a{'^e) — ^ 0(9^) a-e. in 
Qt and with the weak convergence AA((/3e) — t- AA((/3) in L'^{Qt) we are led 
to 

^''((/9e)me(v3£) div77da;dt — > I {ip)m{ip) div rj dx dt and 

Jqt 



\J a{ip£)AA{ipi;)mi;{ipi;) div rj dxdt — > / \/ a{ip) AA{ip)'m{ip) div rj dxdt. 

Jqt 

The next step is to show that m'^{ipe)V^e — ^ m'{ip)Vip in L'^{Qt)- To this 
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end we split the integral in the following way: 

\m'^{(pir)Vipe — m! {ip)V^p\^ dx dt 

Qt 

|m^((/9e)V(/7£ — m'{Lp)V^p\'^ dx dt 

'Tn{|(^|<i} 



+ / |m^(c/?£)V(/?e — m'((/3)V(/9p drc . 

'Tn{|</p|=i} 



Since V^p = a.e. on the set = 1}, see for example Gilbarg and 
Trudinger |GTOH Lem. 7.7], we obtain 



/ \m'^{ips)^'^e — m! {ip)V ip\^ dx dt 

jQTn{\v\=l} 

= I \m!e{'^e)^ Ve\^ dx dt 

<C f \V^e?dxdt^C f iVc^l^drrdt = 0. 

jQti^\\^\=i) JQTn||(p|=i) 



nn{|</3|=i} J QTC^{\^p\=l} 

Although is not continuous, we can conclude on the set {\ipe\ < 1} the 
convergence m'^{ip^) — m'{ip) a.e. in Qt- Indeed, for a point {x,t) G Qt 
with |(/3(x,t)| < 1 and — t- ip{x,t), it holds that < 1 — 6 for 

some 6 > and e small enough and in that region and m' are continuous. 
Hence we have 

m'^{ip£)Vip£ — > m'{ip)Vip a.e. in Qt (3.22) 
and the generalized Lebesgue convergence theorem now gives 

|m^((/7£)V(/3e — m!{ip)Vip\^ dx dt — > , 



/ 

JQ 



Tr\{\f\<i} 



which proves finally m^((^£)V(/9e — )• m'{ip)'Vip in L'^{Qt)- Similarly as above, 
together with the convergences "^'{ipe) — ^'{<f), aife) — a(y') a.e. in Qt 
and with the weak convergence AA^ip^) — t- AA{ip) in L'^{Qt) we are led to 

^'{iPe)m'£{<pe)Vip£ -ridxdt 

^' {ip)m' {ip)V ip ■ T] dx dt and 

't 

^/a{^Pe)AA{ipe)■m'^{'Pe)V<p£ ■ T] dx dt 



I 

jQi 



\J a{ip) AA{ip)'m' {p)^/ ip ■ r} dx dt. 
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Now we are left to show that the second term of the right side in (|3.2ip 
converges to zero. To this end, we spht it in the fohowing way: 

= e "if'i^iipe) d[v{me{(ps)r]) dx dt 

+ e ^^{ipe) dw{me{ipe)v) dx dt 

= ■■ il)e + . 

On the set {\(pe\ < 1 - e} we use that ^[^{(fe) = ln(l + ips) - ln(l -(pe) + 2 
and therefore < | lne| + C to deduce that 

\{I)e\ <e{\ln£\+C) / \diY{me{ipe)v)\dxdt — ^0. 
JQt 

On the set {\fe\ > ^ — we use that me((/9e) = e(2 — e) to deduce 



{II)e=e\2-e) [ ^^iipe)dWridxdt 

<CeVi^i^,)\\LHQT) 

since the last term in brackets is bounded by the energy estimate form 
Lemma 13.71 

For the relation of J and J we note that due to Jg ^ J, ^ J in 
L'^{Qt), Je = v^^el^e) Je and y^mjjp^ — )• m{ip) a.e. in Qt from ()3.22p 
we can conclude 



From the weak convergence ^ J in LP'{Qt) we can conclude that 

/ I (ix dr < liminf / mi.{Lp^)\V dx dr 
jQ(S,t) JQ(S,t) 

for all < s < 4 < T and this is enough to proceed as in Abels, Depner and 
Garcke |ADG12) to show the energy estimate. 

Finally we just remark that the continuity properties and the initial 
conditions can be derived with the same arguments as in [ADG12[ Sec. 5.2, 
5.3], so that altogether we proved Theorem 13.51 



20 



Acknowledgement 

This work was supported by the SPP 1506 "Transport Processes at Fluidic 
Interfaces" of the German Science Foundation (DFG) through the grant GA 
695/6-1. The support is gratefully acknowledged. 



References 



[Abe09a] 



[Abe09b] 



[Abel2] 



[ADG12] 



Abels H., Existence of weak solutions for a diffuse interface 
model for viscous, incompressible fluids with general densities, 
Comm. Math. Phys., vol. 289 (2009), p.45-73. 



[AGG12] 



[ADGK12] 

[Boy99] 

[Boy02] 

[CT94] 
[Che02] 



Abels H., On a diffuse interface model for two-phase flows 
of viscous, incompressible fluids with matched densities, Arch. 
Rat. Mech. Anal., vol. 194 (2009), p. 463-506. 

Abels H., Strong Well-Posedness of a Diffuse Interface Model 
for a Viscous, Quasi-Incompressible Two-Phase Flow, SIAM J. 
Math. Anal, vol. 44, no. 1 (2012), p. 316-340. 

Abels H., Depner D., Garcke H., Existence of weak solutions for 
a diffuse interface model for two-phase flows of incompressible 
fluids with different densities, larXiv:rill.2493t to appear in J. 
Math. Fluid Mech. (2012). 

Abels H., Garcke H., Griin G., Thermodynamically consistent, 
frame indifferent diffuse interface models for incompressible 
two-phase flows with different densities. Math. Models Meth. 
Appl. Sci., vol 22, no.3 (2012). 

Aki G., Dreyer W., Giesselmann J., Kraus C., A quasi- 
incompressible diffuse interface model with phase transition, 
WIAS preprint no. 1726, Berlin (2012). 

Boyer F., Mathematical study of multi-phase flow under shear 
through order parameter formulation, Asymptot. Anal., vol. 20, 
no.2 (1999), p.175-212. 

Boyer F., A theoretical and numerical model for the study of 
incompressible mixture flows, Comput. Fluids, vol. 31 (2002), 
p. 41-68. 

Cahn J. W., Taylor J. E., Surface motion by surface diffusion. 
Acta Metall., vol. 42 (1994), p. 1045-1063. 

Chen L.-Q., Phase-field models for micro structure evolution, 
Annu. Rev. Mater. Res., vol. 32 (2002), p. 113-140. 



21 



[DU77] Diestel J., Uhl Jr. J. J., Vector Measures, Amer. Math. Soc, 
Providence, RI, 1977. 

[DSS07] Ding H., Spelt P.D.M., Shu C, Diffuse interface model for in- 
compressible two-phase flows with large density ratios, J. Comp. 
Phys., vol. 22 (2007), p.2078-2095. 

[EG96] Elliott CM., Garcke H., On the Cahn-Hilliard equation with 
degenerate mobility, SIAM J. Math. Anal., vol. 27, no. 2 (1996), 
p 404-423. 

[Gru95] Griin G., Degenerate parabolic equations of fourth order and a 
plasticity model with nonlocal hardening, Z. Anal. Anwendun- 
gen, vol. 14 (1995), p.541-573. 

[GPV96] Gurtin M.E., Pohgnone D., Vinals J., Two-phase binary fluids 
and immiscible fluids described by an order parameter. Math. 
Models Meth. Appl. Sci., vol. 6, no.6 (1996), p.815-831. 

[GTOl] D. Gilbarg, N. S. Trudinger, Elliptic PaHial Differential Equa- 
tions of Second Order, Springer 2001. 

[Hil70] Hilliard J. E., Spinodal decomposition, in Phase Transforma- 
tions, American Society for Metals, Cleveland, 1970, p.497-560. 

[HH77] Hohenberg P.O., Halperin B.I., Theory of dynamic critical phe- 
nomena. Rev. Mod. Phys., vol. 49 (1977), p.435-479. 

[Lio69] Lions J.-L., Quelques Methodes de Resolution des Problemes 
aux Limites Non lineaires, Dunod, Paris, 1969. 

[LT98] Lowengrub J., Truskinovsky L., Quasi-incompressible Cahn- 

Hilliard fluids and topological transitions, R. Soc. Lond. Proc. 
Ser. A Math. Phys. Eng. Sci., vol. 454 (1998), p.2617-2654. 

[Mue85] Miiller I., Thermodynamics, Pitman Advanced Publishing Pro- 
gram. XVII, Boston-London-Melbourne, 1985. 

[Rou90] Roubicck T., A generalization of the Lions-Temam compact 
embedding theorem, Casopis Pest. Mat., vol. 115, no.4 (1990), 
p.338-342. 

[Sim87] Simon J., Compact sets in the space U'{f),T;B), Ann. Mat. 
Pura Appl. (4), vol. 146 (1987), p.65-96. 

[SohOl] Sohr H., The Navier-Stokes Equations, Birkhauser Advanced 
Texts: Basler Lehrbiicher, Birkhauser Verlag, Basel, 2001. 



22 



